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The particle-particle hole-hole ring-diagram summation method is employed to obtain the equa- 
tion of state of asymmetric nuclear matter over a wide range of asymmetry fraction. Compared 
with Brueckner Hartree-Fock and model-space Brueckner Hartree-Fock calculations, this approach 
gives a softer equation of state, increased symmetry energy and a lower value for the incompress- 
ibility modulus which agrees quite well with the values used in the hydrodynamical model for the 
supernovae explosion. 



1. INTRODUCTION 

A primary aim of microscopic nuclear theories is to derive the various nuclear properties such as the binding energy 
per nucleon (BE/A) and saturation density (po) of nuclear matter, starting from fundamental nucleon-nucleon (NN) 
interactions. The well-known BHF approach is one such standard nuclear matter theory. In terms of G-matrix 
diagrams, the BHF theory is, however, only a lowest-order approximation; the ground-state energy shift AEq for 
nuclear matter, due to the NN interaction, is given merely by the first-order G-matrix diagram, fig. 1(a), namely 

AE* HF = J2^ a n b (ab\G BHF (cj = e a + e b )\ab) (1) 

ab 

In the above the n's are the unperturbed Fermi-Dirac distribution functions, rik=\ if k < Uf and =0 if k > kF- kp 
is the Fermi momentum. 




FIG. 1: Some lower order ring diagrams 

The single-particle (s.p.) energies denoted by e, are determined self-consistently using the BHF theory. As is well 
known, this approach has not been very successful; it has, in general, not been able to simultaneously reproduce the 
binding energy per nucleon (BE/ A = — 16 ± 1 MeV) and the saturation Fermi momentum (fcf? 4 = 1.35 ± 0.05 fm -1 ). 
When plotted on a energy-density plane, results of various BHF calculations for nuclear matter invariably lie, more 
or less, on a band, the Coester bandi, which significantly misses the "experimental (empirical) box" for BE /A and 
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The ground state of nuclear matter in BHF theories is treated as merely a Fermi sea. Particle-hole fluctuations 
near the Fermi sea which represent the long-range correlations are not considered. It should be of interest to allow for 
such Fermi-sea fluctuations, and they may be important in determining the stiffness of the nuclear-matter equation of 
state (EOS). Yang, Heyer and Kuc 2 proposed an elegant and rigorous method for summing up particle-particle hole- 
hole (pphh) and particel-hole (ph) ring diagrams to all orders for the calculation of ground state energies of general 
many-body system. Inclusion of these classes of diagrams to all orders takes into account the Fermi sea fluctuations 
and long-range nuclear correlations. 

With this motivation, several calculations for symmetric nuclear matter have been carried out with the inclusion of 
certain class of ring diagrams to all order o 3 i i 5 . 

In comparison with conventional BHF 6 ' 7 calculations of nuclear matter, the inclusion of the particle-particle hole- 
hole (pphh) ring diagrams to all orders has the desirable effect of both increasing the nuclear matter binding energy 
and lowering its saturation density. The final expression for the ground-state energy shift in the pphh ring diagram 
summation with a model-space reaction G- matrix (G M ) is given as 

AC 9 = f Q ^tv{Y m (X)Y^(X)G M (w = A^- 2 (A))A} (2) 

Comparing with the corresponding BHF result of eq. the main difference between the two methods is the 

replacement of the unperturbed occupation factors n, BHF G-matrix G BHF and starting energy (e a + e&) in the BHF 
expression by the RPA amplitudes Y, model-space G-matrix G M and RPA energies A, respectively. In the above A 
is a strength parameter introduced to facilitate the calculation, as will be discussed later. 

In the present paper, we would like to extend the above ring-diagram scheme to asymmetric nuclear matter, which 
is in several ways of more physical importance than symmetric (N=Z) nuclear matter. The study of the EOS of 
asymmetric nuclear matter has become, in the last few years, a subject of renewed interest particularly in connection 
with astrophysical problems^, such as supernovae explosions and the evolution of neutron stars. For these physical 
processes, the nuclear matter involved is predominantly not symmetric, and it is the EOS of asymmetric nuclear 
matter (with a large neutron excess) which plays an important role for them. Furthermore, the nuclear matter probed 
by heavy-ion experiments is also generally asymmetric. Here the proton- neutron ratio is about 2/3 for both the target 
and the projectile, and thus the resulting nuclear matter is likely to be asymmetric with the same proton-neutron 
ratio. The range of densities sampled by astrophysical systems such as the supernovae and the neutron star vary over 
an appreciably wide range as do their isospin asymmetries. 

According to the model of prompt explosion 10 , which has been widely employed to explain the explosion mechanism 
of a supernova, the nuclear-matter EOS needs to be relatively soft& for the explosion to take place. It should be of 
interest to investigate the effect of ring diagrams on the stiffness of asymmetric nuclear matter, as we shall carry out 
later. To our knowledge, such investigations have not yet been performed. 

Asymmetric nuclear matter calculations have been done using the Skyrme interactionsii*i^, the Gogny interactioni 3 - 
and using the Brueckner-Bethe-Goldstone (BBG) approachi^i^. Our present calculation is a continuation of the 
model-space BHF (MBHF) calculations for asymmetric nuclear matter carried out by Song, Wang and KuoA&. In 
the past, the asymmetric nuclear matter properties were often extracted by interpolating the two extreme cases of 
symmetric and pure neutron matter, with an empirical parabolic approximationiii^. The validity of this empirical 
practice seems to have not been investigated, and we would like to carry out such an investigation by carrying out 
a sequence of ring-diagram nuclear matter calculations, covering a wide range of proton-neutron ratios and baryon 
densities. 



2. FORMALISM 



Asymmetric nuclear matter is a system consisting of N neutrons an Z protons with N ^ Z. For symmetric nuclear 
matter N and Z are identical with the same Fermi momenta. For asymmetric matter however, the neutrons and 
protons are treated as non-identical particles with different Fermi mommenta. 

We introduce a parameter a as a measure of the asymmetry in nuclear matter, namely 

(Pn - P P ) (N-Z) 

a = = — -a — \A) 

P A 

where p, p n and p p are respectively the nuclear, the neutron and the proton densities. The neutron and the proton 
Fermi momentum are 

k n F = (37r 2 p„) 1/3 k F = (3tt 2 p p ) 1/3 (4) 
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An average Fermi momentum is defined as 

o 1/3 

k F = (-n 2 p) (5) 

3. MODEL-SPACE G-MATRIX AND SINGLE-PARTICLE SPECTRUM 

As with the symmetric case we start with a Hamiltonian H = T + V . Introducing a single particle (s. p.) potential 
U we rewrite it as H = (T + U) + (V — U). V is the NN potential such as the Parish or the BonnSi potential. A 
model space P is defined as a configuration space where all nucleons are restricted to have the momentum k < kM, 
ku being the momentum space boundary of P. A typical value for /cm is 2k F where k F is the Fermi momentum. As 
we shall discuss later, most of the calculations in the present work are performed using /cj\/=3.0 fin . 

Similar to the case of symmetric nuclear matter we use a model-space Hartree-Fock methodSLSSiS 3 . to determine U. 
This leads to the following self-consistent equations for the s.p. spectrum ef' 1 , namely 

e f - t ki +r ki {h, n ) (6) 

T ki {cj, n) = 2 J2 (hkj\G M \w + ej)\kikj) h < k M (7) 

r fei (w,ri) - h>k M (8) 

Note that the subscript i represents both momentum and isospin, namely i = (kj, r^) with kj denoting the momentum 
and Ti the z component of the isospin of the ith nucleon. The single particle kinetic energy is t ki — h 2 kf /2m, G M is 
the reaction matrix to be specified later. The Fermi momentum is represented by kp with neutron Fermi momentum 
k F for Tj — — o, and proton Fermi momentum k F for r,j — ^. The model-space momentum boundary is fcju and is 
taken to be the same for neutrons and protons. The s.p. potential is the one-body vertex function T evaluated at the 
self-consistent energy u> = e ki 

U(ki,Ti) =T ki (eh i ,T i ) , i = n,p (9) 

U{ki,Ti) is determined self-consistently for ki < k^.j, and for ki > ku we set U = 0. We also use an effective mass 
description for the single particle spectrum as 



M = / (h 2 /2m* q )kf + A 9 , h < k M 
(h 2 /2m)k 2 , k > k M 



with four parameters m*,A qi (q = n,p). The effective mass m* and the zero point energy A are determined self- 
consistently. 

The model-space reaction matrix G M satisfies the Bethe-Goldstone equation 
/ • • i siM i \ i v \ (ij\V\kl)Q M (k,l)(kl\G M {u;)\mn) 

{i]\G M (u:)\mn) = (ij\V\mn) + > — — — — — — ! '- (11) 

*—f w-e k -ei 

k,l 

In the above equation i, j, m, n, k and I are single particle states, each with momentum k, isospin r. V is the nucleon- 
nucleon interaction. The energy variable u> in the denominator of eq. (11) is given by 

a, = ef + ef (12) 

The two particle correlations considered are those where at least one particle is out of the model space. Hence the 
Pauli operator Q M is defined as 

nM (u i\ — / !> if max(fc fc , fc ; ) > k M and k k > k T F k ,ki > hp 1 , . 

Q ^ fe ' L > - \ 0, otherwise 

Note that Q M is different for each of the following three cases of the intermediate two nucleon state. 

For the nn (pp) case (figEJa)), the intermediate nucleons are identical, and only k^[ and k q F , q — n or p, enter 
the calculation. For the np or pn case however, kM, k F and k p F all play a role in determining Q M (fig^b)). It is 
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FIG. 2: Intermediate two nucleon states in the asymmetric case (a) nn or pp (b) np or pn. 
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FIG. 3: Pauli operator Q M (k k ,ki 



convenient to carry out the above calculation in the relative and center of mass (RCM) frame. We choose our relative 
momentum k and center of mass momentum K as 



1 



(k k - k,) 



K = (k k + kO 



(14) 



First we replace the Pauli exclusion operator Q M , which is a function of the laboratory momenta by its angle 
average approximation Q M . We divide the plane of the two laboratory momenta into three regions A, B and C as 
shown in fig|31 The values of Q in the three regions are shown in the figure. Each of the regions is transformed into 
the RCM frame. Then the angle-averaged value of Q M is 



Q 



Of 



Qc 



(15) 



where, for region B, we have 



Q%(k,K) = 





((k 



\K) ( _ 
{2k 2 + \K 2 - 



k 2 M )/2kK 



k 2 



)/2kK 



{k n F Y)/2kK 



(K - (k n F ) 2 )/2kK 



regions a, 
region c 
region d 
region e 
region / 



(16) 



The subdomains a to e are shown in figQJand 2(fcg) 2 = [(k F ) 2 + k 2 M ]. Angle average approximations are standard 
(and indispensible) in treating Pauli exclusion operators in nuclera matter calculations and are generally considered to 
be accurat e) 24 ! 25 . This technique has the advantage of making the model-space reaction matrix diagonal in the RCM 
vriables. Recent studies of an exact pauli operator calculation have been reported^. It is not very clear whether such 
calculations show an appreciable difference in the binding energy per nucleon and the saturation density from previous 
studies. Angle averages for region C may be obtained from the above by substituting k p F for k F . The calculation 
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(a) (b) 



FIG. 4: Angle averaged Pauli operator and Qa ■ 




above was an illustration of the case where the intermediate state contains one neutron and one proton. The other 
two possibilities may be readily obtained from above by suitable substitution of the relevant Fermi momentum. 
A similar RCM mapping of region A in fig|3| gives 

region a 

region b (17) 
jK 2 — k 2 M )/kK region c 

For symmetric nuclear matter k F = k F so that Qg — and the angle averaged Q M is the same as that given in 

Using the angle averaged Pauli operator the model space reaction matrix can be decomposed into separate partial 
wave channels as 

wg^kvwv) = (kiwwn + 2 - fV ^ 2E (fc^i^^^ M (fe^^)(fe^^ M (^^)i^o (18) 

k Jo yf u-H (k",K) 

where w denotes the partial wave quantum numbers (ll'SJT), and K is the center of mass momentum. The K and 
(SJT) quantum numbers associated with the bra and ket vectors have been suppressed for simplicity. For example, 
(kl\ should in fact be {klSJT,K\. Our convention for plane waves is 

(r\klSJ) = 3l (kr)y lS j(r) (19) 

where ji(kr) is the spherical Besscl function, and y is the vector spherical harmonics corresponding to 1=S+J. The 
angle averaged reaction matrix G M (u>, Kw) is diagonal in K and w. This is a consequence of using angle averages for 
the projection operator Q M and the energy denominator. 

From eq. (|12fl the energy variable in the laboratory frame is e, + e%. Using RCM variables, the energy denominator 
u> for the neutron spectrum is given by 



& , 2 & , A 

= — : k + + A « + a p " 

m* 4m* 



2m* 2m* 



(k n ) 2 (20) 



Similar expression is obtained for the proton spectrum by replacing the subscripts and superscript n by p. The 
other term in the denominator Hq = e m + e„ is the unperturbed energy of the intermediate states and is also angle 
dependent. The momentum variables k m and k n corresponding to k and K may be in either of the three regions A, 
B or C. In region A, both intermediate particles have momentum larger than ku- Therefore we have 

H £(k,K) = -k 2 + ^K 2 (21) 
u m 4m 

In region B, we have a proton with momentum larger than ku and a neutron with momentum between k F and fcjvf- 
We use an angle average approximation for k 2 ^ i.e., (fc^) = \ [{k F ) 2 + k 2 M ] and obtain 

H*(k, K) = h -k 2 + ^K 2 + An + f[^- l][{k F ) 2 + k 2 M ] (22) 
m 4m 4m m* 



G 



The spectrum of Hq for region C is the same as that of region B with the subscripts and superscript n changed to p. 
Finally the single particle potential in the (RCM) frame for ki < kp with (q = n 7 p) is given as 



U(ki,r q ) = - £ ( 2J + X ) / dk k2G isjr qT] (K Ki) 



ISJtj 

~ ( 2J + 1) dk fc t(fc?) 2 - + ik(ki - k)]Gf sjTqrj {k,K 2 ) 

ISJTj 



k- 



(23) 



And for k\ > kp 1 but less than ku we have 

*i> /+i) £ 



ISJTj 



tf(*i,T,) = 4" E ( 2J + !) / * fc t( fc ?) 2 - *? +4fc(fci - k)]G$ Jr<T; ,(fc,tf 2 ) (24) 



where 

fc+ = i(fcp +fcl) 

- 4(^ + fc 2 ) 

Kl = 4(kf + k 2 )-(2k + k 1 -k T F !)(2k + k 1 + k^) (25) 
and the partial wave reaction matrix elements are given by 

Gf§ jTiT .(k,K) = {klT l T 1 \G M (u;,KlSJ)\klT 2 T J ) (26) 

Eq. I|18l) is in isospin representation with well defined total isospin T of the two nucleon state. The reaction matrix 
elements in the neutron-proton representation i.e, eq. (|26|) . are related to those in the isospin representation by 
Clebsch- Gordon coefficients; which explains the factor of 2 in the expression for the s.p. potential. 

The s.p. potentials U(ki, r n ), U(ki, t p ) and the spectra e„, e p are calculated self-consistently as described previously. 
A main purpose here is to convert the strong V interaction to a well-behaved G-matrix interaction. In so doing, 
one must make sure that there is no double counting. Thus a double-partitioned approach is adopted, treating the 
nucleon- nucleon correlations with high momentum (i.e. those with Q M = 1) within the G M matrix, while taking care 
of the low-momentum correlations by the pphh ring diagrams. We would like to express the energy shift in terms of 
the model space G-matris*. 

Proceeding as in Ref. (3) we define a model-space G-matrix and formulate an expression for the energy shift in 
terms of the transition amplitudes as 

A£T = / y E E y «^'< X ) Y ^ = ^r 2 H])A (27) 



i>j,k>l 

(A-2) |p 



with 



Ei<W ' + £ i) + (1 - «i - %-)A£ye/M}r ro (e/, A) = A^- 2 (u)Y m (ij, A) 

e>f 



.A-2 



(u) (28) 



For the neutron-neutron (nn) or the proton-proton (pp) cases, the inequality signs for the summation indices i, j, k 7 I 
restrict the momentum ki > kj, similarly for k^ and ki to avoid double counting. The neutron-proton (np) or the 
proton-neutron (pn) case is more complicated. Here a free summation over the indices gives four terms with identical 
contribution (two for each of np and pn cases) and a factor of (1/2) would be needed for correct counting. Thus in our 
calculation for this case we have confined the indices i, fc, e to neutrons and j, Z, / to protons with their momentum 
summations unrestricted. 
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FIG. 5: Neutron (U n ) and proton (Up) s.p. potentials. 
I 




FIG. 6: Single particle spectra for neutron and proton calculated with the model space approach at a = 0.2 and 0.8 for 
&f = 1-35 fm" 1 with Paris interaction. 



4. RESULTS OF THE S.P. SPECTRUM CALCULATION 



Fig.J^J shows the neutron (upper part) and the proton (lower part) s.p. potentials for different a. As the proton 
fraction decreases the s.p. potentials become less deep and U p vanishes when the proton fraction becomes zero. 

For the asymmetric case as well, the spectrum is continuous from momentum to as shown in figEl unlike 
the usual BHF spectrum which has a large discontinuity at fc^- Beyond ku we use a free s.p. spectrum since our 
method is designed for determination of the s.p. spectrum within the model space only. Our spectrum has a small 
discontinuity of around 4 — 5 MeV at ku- 

Table 1 lists the m/m* values and zero point energies at a = 0.2 and a = 0.8 for various average Fermi momentum 
at a model space boundary of 3.0 fm" 1 . With increase in a, A„ increases and A p decreases for small kp. For the 
same asymmetry fraction, both A n and A p become deeper with increase in Uf- The change in the effective mass with 
the asymmetry fraction is very small though for the same asymmetry fraction it increases with kp. For a = 0.0, the 
s.p potential, the effective masses and the zero point energies match the ones calculated previously for the symmetric 
case which also serves as a test for the reliability of our asymmetric matter calculation. The calculated binding energy 
for a given combination of the asymmetry parameter a and the Fermi momentum kp depends on the model-space 
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TABLE 1: Typical m* and A values calculated with the Paris NN interaction at a = 0.2 and 0.8. 



k F (fen -1 ) 


rn/m* n 


A„ (MeV) 


m/rrip 


A p (MeV) 


a = 0.2 1.00 


1.10 


-29.44 


1.10 


-31.45 


1.20 


1.15 


-46.61 


1.16 


-48.89 


1.35 


1.21 


-61.69 


1.21 


-63.73 


1.50 


1.26 


-77.38 


1.26 


-76.83 


1.70 


1.34 


-96.75 


1.34 


-97.13 



a = 0.8 



1.00 


1.08 


-24.66 


1.10 


-33.05 


1.20 


1.13 


-40.05 


1.15 


-49.56 


1.35 


1.17 


-52.16 


1.19 


-60.60 


1.50 


1.23 


-66.47 


1.24 


-74.05 


1.70 


1.32 


-86.55 


1.32 


-91.68 



size. On minimizing the energy against the model-space momentum, Song, Wang and Kuoi& found for their MBHF 
calculations that for each combination of a and kp there was one value of ku for which the energy was a minimum. 
For the same a, the energy minimum shifted towards smaller values of ku- For values of kF ranging between 0.50 to 
1.80 the minimum values of ku were between 2.80 to 3.4. Based on this we have made our choice of ku to be 3.0 
fin- 1 . 



5. THE RPA EQUATION 



Let us recall our RPA-type secular equation 
y"]{(ej + e 3 )Si 3 . e f + {niUj - ninj)\L i: j e f(uj)}Y m (ef, A) = (i m (u, X)Y m (ij, A) 

e>f 

with the self consistent condition 

Hm(u, A) = uj = u> m (\) 
L ijef (uj) ee G™ ef (w) 

The above RPA-type equation is in laboratory momentum variables. As for 
above equation into its RCM representation. We also do an angle average for 



- (n< 



We define a function Qniki, kj) = 1 — (n$ 



It is equal to +1 



1 

kj. Qr is +1 in regions A and B (i,j both are particles); is —1 in region C (i 
for all other regions. The value of Qr depends on the angle between k and K 
of the total isospin T of the two nucleon state i.e., t z = Ti + Tj. As the two 
np (pn), t z can take the values -1,1 or 0. For the case when t z = +1 or — 1 
situation is no different from the symmetric case {Qr is the same as eq. (4.9) 
angle averaged Qr as 



(29) 



(30) 

the symmetric case, we transform the 
the occupation factor (hifbj — riiUj) = 
or - 1 depending on the values of ki and 
j are both holes) and is equal to zero 
Let us denote by r z the ^-component 
nucleon state can be either nn, pp or 
the two nucleons are identical and the 
in Ref. (3). For r 2 = we obtain the 



-1 



Q R {k,K) = I 



-\xi\ 
1 



region 1 
2 
3 

xi\ 4 
x 2 \ 5 
, min(|:ci|, |x 2 |) 6 



(31) 



where 



fc" 4- kl 
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The above angle averages are obtained under the assumption that all values for the angle between k and K are equally 
likely. The regions 1 to 6 refer to the regions in the (k,K) plane as shown in the figure. The replacement of the 





(a) 



2k" 2k M K 



(b) 



FIG. 7: Angle average of occupation factor Qn(ki,kj 



occupation factor Qji by its angle-averaged quantity greatly simplifies the RPA-type secular equation we started with. 
Now it can be decomposed into separate partial wave channels 



f 2k' 2 
J2 / dk'{e kK S(k - k')S w + A Qi 



(k,K)(kl\L(uj,K)\k'l')}Y m (k'l'k,X) 



(uj,X)Y m (klK, A) 



(32) 



where Qa(k,K) is actually Qu(k,K,kp ,kp , few)- is the unperturbed energy. There is a subtle point. For the 

nn or pp cases where the two nucleons are identical this unperturbed energy may be taken simply as i^k 2 + -j^K 2 + 

2U(k\,T q ), where k\ = (k 2 + \K 2 ) is the angle averaged value for the momenta of the two nucleons and q = n or 
p. For the np (pn) case, this averaging is more complicated. One way around this is in the choice of relative k mesh 
points. For each K mesh point we choose the k such that the RCM values of k n and k p are either both in the hole 
region or both in the particle region. With such a choice of mesh points the sum of the squares of the individual 
neutron and proton momenta in the RCM frame, both being either holes or particles is k^ + k 2 — (k 2 + jK 2 ). Hence 
the value of the unperturbed energy is 

h 2 h 2 

—k 2 + —K 2 + U{k x ,T n ) + U{k u r p ) 
m Am 

with k\ 2 = (k 2 + \K 2 )/2. 

The wave function (kl) denotes the RCM partial wave function (MSJT1T2, K) and similarly for (k'V). The above 
equation is to be solved together with the self consistent condition fJ, m (u>, A) = w~(A) giving the self-consistent solution 

The vertex function L(u>, T1T2, K) in the above equation is the irreducible vertex function which has both one-body 
and two-body G-matrix diagrams^. These contribute significantly to the depletion of s. p. orbits below kjp , especially 
at high density. As we are working in the RCM frame an angle average approximation is employed to obtain L in the 
RCM frame for the t z = case as 



(kl\L(u,T n T p ,K)\k'l' 



{klT n T p \G M {uJ,Kw)\k'l'T n T p ) 

+5 kk i8w{Tj, ri (u) - ej. p ,T n ) - U(k n ,r n ) 



(33) 
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where 



k 2 p = (fc 2 + -K 2 )I2 



% = TT~ + U(k p ,T p ) 
p 2m 



The other two cases for t z = — 1 or 1 can be obtained by suitable substitution of t p and r n . 



6. RESULTS AND DISCUSSIONS 



(34) 



6.1. Binding energy 



The energy shift is given as 
G 



W T; 

k M rk 



r 1 r 2kT F 

tEE( 2J+1 ) dX dKK 

\ „, Jo Jo 



AE$ P _ 

A ~ ^ 2 [(^)3 + (^)3] 

x V / dkk 2 I dk' k ' Y£(MK,\) x (klriT2\G M (uj , Kw)\k'l l TiT2)Y m (k'l / K, A) 

where t z — —1, 0, 1 for T = 1 and t z = for T = and 

f kp for r z = -1 

k F = \ (k F + k p F )/2 for t z = 
[ fej for r 2 = 1 

The binding energy per nucleon is given as 

-BE ih\ [{i + am)2 + a)(fc p )2] + A£T 



(35) 



(36) 



(37) 



A 20m "~ ' ~*' v "" ' v ~ ~' v FJ 1 ' A 
We present the results of our ring calculation in table 2 and in figs and fusing Paris NN interactions, for various 




FIG. 8: BE I A from Ring diagram summation at three asymmetry values : a — 0.0, 0.6 and 1.0 
values of the Fermi momentum and the asymmetry fraction. The binding energy at a = 0.0 (symmetric nuclear 
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FIG. 9: BE/A versus Uf with Ring summation and MBHF calculation at a = 0.8. The ring curve shows a saturation which 
is not present in the MBHF curve. 



TABLE 2: BE/A in MeV for various combinations of a and Uf (fm 1 ) with ku = 3.0 (fm 1 ) and Paris interaction. Results 
from (MBHF) calculation are also given for comparison. 



a/k F 0.65 1.00 1.20 1.35 1.50 1.70 



0.0 


-6.26 


-10.89 


-14.02 


-15.68 


-15.66 


-12.76 




(-3.74) 


(-8.60) 


(-11.40) (-13.31) 


(-14.43) 


(-13.47) 


0.2 


-5.87 


-10.49 


-13.47 


-14.97 


-15.00 


-11.64 




(-3.42) 


(-7.94) 


(-10.56) (-12.38) 


(-13.27) 


(-12.12) 


0.4 


-5.78 


-9.60 


-12.32 


-13.30 


-12.76 


-8.58 




(-2.44) 


(-5.99) 


(-8.00) 


(-9.32) 


(-9.60) 


(-7.76) 


0.6 


-5.18 


-7.42 


-9.62 


-10.04 


-8.87 


-2.40 




(-0.91) 


(-2.68) 


(-3.65) 


(-4.02) 


(-3.30) 


(-0.03) 


0.8 


-1.90 


-3.52 


-3.83 


-3.54 


-1.11 


6.40 




(1.29) 


(2.00) 


(2.51) 


(3.52) 


(5.70) 


(11.00) 


1.0 


4.03 


8.11 


11.03 


14.22 


19.02 


29.65 




(4.08) 


(8.07) 


(10.84) 


(13.74) 


(17.85) 


(26.20) 



matter) by our ring calculation is -15.93 MeV which is good agreement with the empirical value. The saturation 
Fermi momentum is 1.43 fm -1 . As the asymmetry fraction increases ( the proton fraction decreases), both the 
binding energy and the saturation Fermi momentum drop till at a = 1.0 (zero proton fraction) there is no saturation. 
An interesting result of the ring calculation is the existence of saturation at a = 0.8 which is not present in the MBHF 
calculation. 
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TABLE 3: E*™ 9 calculated from eq. (38) and eq. (39) and E, 



MBHF 
sym 



kF(fm 


1\ pRing 
) sym 


pMBHF 
^ sym 


Empirical 


0.65 


9.84 


8.17 


10.29 


1.00 


18.48 


16.07 


19.00 


1.20 


23.81 


21.23 


25.05 


1.35 


28.29 


25.42 


29.90 


1.50 


32.79 


29.72 


34.68 


1.70 


40.53 


35.61 


42.41 



6.2. Symmetry energy 



The symmetry energy is defined as 
1 d 2 W(a,k F ), 



da 2 



(38) 



where W(a, k F ) is the binding energy at a given a and k F . To ascertain the nature of the dependence of the binding 
energy on the asymmtery parameter we have tried fitting a polynomial curve of leading order a 2 and higher for the 
residuals upto to the order of a 6 to our data. Using eq. I|38l and the parameters of the above fit we obtain the E£££f 
for various values of k F as shown in table 3. The corresponding values obtained by the MBHF calculation 1 ^ are also 
shown for comparison. The value of the symmetry energy calculated by our method is consitently higher than the 
corresponding MBHF calculations as is shown by fig llUI Thus the ring diagram summation improves the symmetry 
energy values as well. At the saturation density(po = 0.17 fm~ 3 ) E^™ 9 ls 28 MeV which is close to the value of 
31 Mev reported in an independent calculation^. Empirically^,^™ 9 may be simply evaluated from the two extreme 
cases of pure neutron matter and symmetric nuclear matter shown in the fourth column of table 3. 



E sym (k F ) = W(l,k F ) -W{0,k F ) 

The E sym calculated with eq. (2EJl and eq. (2H|) differ by about .5 MeV to 1.8 MeV . 



(39) 



35 

> 

<D 

5 

E 30 - 

I 

LU 

g? 25 



'esymjmbhf.daf 
'esym_ring.dat' 



1 1.2 1.4 

Fermi momentum k_F 



FIG. 10: Symmetry energy obtained from Ring (upper curve) and MBHF (lower curve) calculations. 
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TABLE 4: kq(o) calculated from eq. (40). 
a k s F a \fm- L ) ttpja) Wo(a) 



0.0 


1.43 


114.51 


-15.93 


0.2 


1.43 


100.20 


-15.25 


0.4 


1.39 


92.74 


-13.33 


0.6 


1.33 


67.35 


-10.05 


0.8 


1.24 


22.95 


-3.85 



6.3. Incompressibility 

In describing the iron-core collapse of a presupernova using hydrodynamical models, the basic physical inputs are 
the initial mass of the core and the nuclear EOS. The EOS of asymmetric nuclear matter exhibits a minimum which 
disappears before pure neutron matter is reached. Therefore we expect that the incompressibility modulus decreases 
with respect to the symmetric case and vanishes before the proton fraction vanishes. The incompressibility modulus 
is defined as 

/ \ 72 i \ d 2 W(a,k F ) 

«o(«) = 4 («) " \k F =k Fo ( a) (40) 

where (&) is the saturation Fermi momentum at the given a. One of the most sophisticated investigation of the a 
dependence of ko and po so far has been done in Ref. (14) in the framework of BBG (Brueckner-Bethe- Goldstone) 
theory They found that for low a values (a < 0.4) kq and po show a parabolic dependence on a 

K {a) = Ko (0)(l-aa 2 ) (41) 
p (a) = p (0)(l-fo« 2 ) (42) 

with k (0) = 185 MeV, a = 2.027 and p Q (0) = 0.289 fm -3 , b = 1.115. In table 4 we give the k values obtained 
from the ring calculation. Our k values also seem to obey a parabolic dependence on a but with a = 1.21 and 
kq(0) — 112 MeV. We may note here that the values of Ko(0) and a in Ref. (14) were obtained by means of a least 
square polynomial fit to the BHF values of the binding energy and the values are quite sensitive to the degree of the 
polynomial used for the fit. 

Our value of kq(0) = 112 MeV is in good agreement with the one suggested by Brown and Osnes 2 ^ for symmetric 
nuclear matter. At a = 0.33 the value of kq from our fit is 97 MeV which is close to the empirical value of 90 MeV 
used by Baron, Cooperstein and Kahana&to get the maximum explosion energy in their hydrodynamical calculations. 
From their numerical investigation the authors of Ref. (8) concluded that the softening of the EOS plays a crucial 
role in generating prompt explosion for stars in the mass range of 12 — 15 M Q (where Mq ~ 2 x 10 23 g is the sun's 
mass). This conclusion has been questioned by more refined calculations^iffi^i. However, even in models where the 
direct explosion mechanism fails, a softer EOS is helpful to the shock. The decrease in incompressibility with increase 
in a is quite intuitive when we consider that going from bound symmetric matter to unbound neutron matter, the 
minimum of the binding energy gradually disappears. Therefore from eq. (40) the nuclear incompressibility decreases 
and vanishes for a certain value of a. 



Conclusion 



In conclusion, a fully microscopic calculation has been done using the ring-diagram summation for the EOS of 
asymmetric nuclear matter. The numerical computation has been done with both the Parish and the BonnSl po- 
tentials and the results are in satisfactory agreement with the empirical values. The model-space size is treated as 
a parameter. The inclusion of ring-diagrams increases the role of tensor forces and both the binding energy and the 
saturation density values are lower and closer to the empirical values than those obtained with previous calculations. 
The symmetry energy values also show an improvement, and is in better agreement with other independent cal- 
culations. The derived incompressibility exhibits a parabolic dependence on the asymmetry fraction, in qualitative 
agreement with the empirical asymmetry dependence used in the literature. 

The ring diagram approach employs an infinite order summation of particle-particle hole-hole ring diagrams. The 
infinite oder summation technique of Yang,Heyer and Kuo 2 " is applicable to particle-hole ring diagrams as well. But 
these diagrams have not been included in nuclear matter calculations. There is reason to believe that the effect of 
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the particle-hole ring diagrams on binding energy calculatins is not very appreciable and they are less important than 
the particle-particle ring diagrams. The lowest-order ring diagram is the pphh diagram of first order in G M . The 
second-order diagram may be taken as either pphh or ph and is second order in G M . Thus the contribution to the 
ground state energy shift from the ph diagrams comes from the three vertex diagram which is third order in G M (fig. 
1(d)). Studie o 32 i 33 i 34 have indicated that the particle-hole diagrams converge rapidly and they may not be impotant 
in nuclear matterbinding energy calculations. This view is also supported by some Lipkin model calculations^. 

As noted earlier, an interesting result of our ring calculation is the existence of saturation at a — 0.8 which is not 
present in the MBHF calculation. This behaviour would be of relevance in astrophysical systems which are essentially 
neutron rich. It would interesting to study the neutron star properties with this method. To do that one needs to 
extend the present calculation to higher densities in a relativistic frame work. 

The authors would like to thank Prof. G. E. Brown for his support and encouragement throughout the course of 
this work. 
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